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Abstract 

Traditionally, the diffraction of a scalar wave satisfying Helmholtz 
equation through an aperture on an otherwise black screen can be 
solved approximately by Kirchhoff's integral over the aperture. Ru- 
binowicz, on the other hand, was able to split the solution into two 
parts: one is the geometrical part that appears only in the geometrical 
illuminated region, and the other representing the reflected wave is a 
line-integral along the edge of the aperture. However, this decompo- 
sition is not entirely satisfactory in the sense that the two separated 
fields are discontinuous at the boundary of the illuminated region. 
Also, the functional form of the line-integral is not what one would 
expect an ordinary reflection wave should be due to some confusing 
factors in the integrand. Finally, the boundary conditions on the 
screen imposed by Kirchhoff's approximation are mathematically in- 
consistent, and therefore, rigorously, this decomposition formulation 
must be slightly modified by taking into account the correct B.C.s. 
In this thesis, we use the consistent boundary conditions to derive 
a slightly different decomposition formula which shows that the be- 
havior of the diffracted wave at the edge is exactly just like an or- 
dinary reflection — realizing the conjecture of Thomas Young in the 
18 t/l century. We also derived another decomposition formula which 
avoids mathematical discontinuity encountered by Rubinowicz. In the 
last section we demonstrate that our solution is consistent with that 
obtained by Sommerfeld in the rigorous 2-D plane-wave diffraction 
problem, so our formulation in this sense may describe more accu- 
rately the behavior of diffracted wave near the edge of the aperture 
than Kirchhoff's formula. 
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1 Introduction 



1.1 Diffraction Integral Formulae 

The diffraction of a scalar wave if> satisfying the Helmholtz equation 

(V 2 + k 2 ) i> = 
has the solution (at the field point ff): 

where the integral is performed on a closed surface S which does not enclose 
any source of if), and h is the outward normal of S. Here, G is the Green's 
function satisfying 

(V 2 + k 2 ) G = -<5 (3) (f-r». 

For the problem of diffraction through an aperture on an infinite screen, 
one usually defines S to be the union of the aperture, the screen, and the 
infinity. If the aperture is finite and the screen is opaque, one expects that if) 
decreases fast to zero at infinity, and thus one only has to evaluate the integral 
on the aperture and the screen. Kirchhoff further assumed the boundary 
conditions 

dip dip s 

ip = ips and — — = — — on the aperture, 

4 n dn ( 2 ) 

if) = and — — = on the screen, 
on 

where ip s is the unperturbed source field. Also, he Kirchhoff the Green's 
function 



G = Gk = 



1 e ik ¥- r M 



An \\r — ff\\ 

With these assumptions, Equation ([TJ can be reduced to 

However, the boundary conditions imposed by Kirchhoff is mathemat- 
ically inconsistent, thought it gives good approximations near the axis at 
far field zone. Sommerfeld, on the other hand, suggested another consistent 
boundary conditions 

if) — if> s , on the aperture, , , 

if) = 0, on the screen, 
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and adopted Green's function of Dirichlet type which vanishes on the aper- 
ture and the screen. For example, for a planar screen, 

1 ik 1 1 r~r f || -I p^'H^-^ll 

/~f /~i /~i 

47r ||r — r/|| 47r ||r — fj\\ 

where fj- is the mirror image of ff with respect to the screen. With these 
modifications, Equation (pQ) is reduced to 

if) (f f ) = - j ip s ^-^-da. (5) 

J aperture "'l 

1.2 Maggi-Rubinowicz' Decomposition 

With the Kirchhoff integral formula (Equation fl3])), Rubinowicz was able to 
decompose the field ip (ff) into two parts [I]: one that appears only in the 
ordinary geometrical illuminated region is the source field evaluated at the 
field point f/; the other one is a line integral along the edge of the aperture: 

i>s (ff) ,ff £ illuminated region If e lkp f ( p s x pf \ 

otherwise 47r/ edge P/ Vl + Ps-P// 

(6) 

This formula applies to two special cases. The first one is the plane-wave- 
incidence case in which the source field is a plane wave: ip s (f) = e lkp , where 
p is the distance measured from f to a constant phase plane of the incident 
field. If f happens to lie on the edge of the aperture, then we denote p by 
p s . Also, we denote p s to be the unit vector in the direction of propagation 
of the incident field, as shown in Figure (1). The illuminated region is an 
oblique cylinder, as predicted by geometrical optics. Inside the line integral, 
Pf = \\f— ff\\ is the distance from ff to the edge f, and pf is the unit vector 
of r*— ff. 

The second one is called the point-source-incidence case, in which the 
source field is a spherical wave: ip s = e tkp /p, where p = \\f—f s \\ is the 
distance measured from f* to the position of the point source f s . As before, if 
f happens to lie on the edge of the aperture, then we denote p by p s , and use 
p s to denote the unit vector of f — f s . As expected, the illuminated region in 
this case is an oblique cone with vertex at f s , as shown in Figure (2). 
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Image 




Figure (1) Figure (2) 



This decomposition realized Young's interpretation for the diffraction 
phenomena: Young had once proposed that part of the incident field, which 
is called the reflected field, undergone a kind of reflection at the edge of the 
aperture, and the rest of the incident field, called geometrical field, just went 
through the aperture without any perturbation, and, the final diffraction 
wave was the interference of the two fields [2J. 

But as mentioned earlier, Equation ([3D was derived based on inconsistent 
boundary conditions. Now, with the B.C.s proposed by Sommerfeld (Equa- 
tion (J4])), we are able to modify Equation (J6j) and obtain another similar 
expression which is not only mathematically-consistent, but also preserves 
Young's "field-splitting" interpretation. 

Furthermore, Equation OH]) is not entirely satisfactory in the sense that the 
two separated fields are not continuous at the boundary of the illuminated 
region. This problem had been discussed by John. S Asvestas, and he also 
gave an elegant "solid-angle" representation of ip which successfully avoided 
such discontinuity. However, Asvestas' work was also based on Kirchhoff's 
formula, and therefore the mathematical inconsistency still exists. Besides 
that, the solid- angle representation derived by Asvestas does not reduce to 
electrostatic case in the long wavelength limit k — > 0. In this paper, we 
will give a slightly different decomposition formula with consistent boundary 
conditions and generalize the electrostatic result to the diffraction problem. 

Another unsatisfactory feature of Rubinowicz' decomposition formula 
Equation OH]) is that the functional form of the line integral is not what 
one would expect an ordinary reflection wave should be, due to some con- 
fusing factors in the integrand. And therefore in this paper, we will also 
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derive a neater representation of the line integral that mimics the behavior 
of ordinary reflection in geometrical optics. 

In the last section, we'll explain why the boundary conditions (Equation 
P|) is more suitable by compare the result of our work with Sommerfeld's 2-D 
straight-edge diffraction problem, which is the one of the few examples where 
the Helmholtz equation has an exact solution. We'll see that Kirchhoff's 
solution has a relative deviation from the exact solution. 



2 Modified Expression for Rubinowicz' De- 
composition Formula 

With the boundary conditions Equation we begin from Sommerfeld's 
integral formula Equation (JSJ). Since Gd (dip s /dn) = on the aperture, we 
add it back into the integrand and rearrange a little: 

1> (ff) = - [ (V*^f - G D ^A da 

J aperture V ulb ulb J 

( dG K dtp s \ , f ( i dG* K d^ s \ 
I i) s — K ~K7 ) da + / ( Gk ~fT da > 

aperture V C>n C>n J -/aperture V ^ n C ' n J 

(7) 

and in the last line we identify the first integral is nothing but Equation (j^J) 
and thus equal to Equation The second integral has exactly the same 
functional form as the first one, except for the replacement ff — > r*. So the 
final result is 



ip s (ff) , ff G illuminated region 
, otherwise 



Li _ P.*r' { ,,,, 

47T/edge ' Pf \1 + Ps ' Pf 1 + Ps ' Pf 



(Note that is on the opposite side of the screen, so it always lies outside of 
the illuminated region, and thus there is no corresponding geometrical field.) 
Here pj. is the unit vector of f — ff. The image term G* K contributes another 
line integral to the final expression, which now seems more ugly. Actually, 
Rubinowicz had also derived this result in his paper in 1917. However, for 
some reason he seemed to abandon this result and used Equation <^ in 
his successive papers. In the following sections, we'll show that, with some 
deformation, Equation (jHJ) (or equivalently, Equation (j^J)) can take another 
form which has some merits mentioned in the introduction. 
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3 "Reflective" Representation 



3.1 Motivations from Geometrical Optics 

In ordinary geometrical optics, the reflection phenomena can be compre- 
hended this way: given a source distribution, one draw the "image source" 
behind the "mirror," as shown in Figure (3), where there is a point real 
source labeled by S, and the reflected field is equal to the incident field from 
the image point S*. 




So, if what Young really meant (in the early days when he saw the diffrac- 
tion phenomena) by "reflection at the edge" was the reflection in geometrical 
optics, then we expect the line integral should take the form 

Reflection ~ f e lk ( ps+p f) for plane wave, (9) 

J edge 

^reflection ~ f : for point source. (10) 

J edge Ps + Pf 

Namely, as shown in Figure (4), we imagine that the screen has a finite 
thickness, and as the incident field ip s comes in, it is reflected by the "cut" 
around the aperture, and thus when the reflected field reaches the field point 
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ff, the incident field ip a has propagated for a total optical length p s + pf, and 
therefore ^reflection should take the form as Equation (jUj) or ffTUj) . 

These expectations actually can be accomplished by some deformations 
of Equation dHJ), but let's do it another way: to derive the ^reflection from the 
beginning Equation (J5J), and this will make derivation more neater. 

3.2 Reflection at the Boundary 

Simplify Equation (J5]) a little, and we get 



where r = ||r — r/|| , and let e z be the inward unit normal to the screen 
(e 2 = — n). To simplify the result, we consider the following two cases: 
Case 1. Plane Wave Diffraction 



We assume the wave propagates in the direction perpendicular to the 
screen-and it is reasonable to make this assumption since experimentally it is 
the most common configuration. Under this postulation, ip s = e tkps where p s 
is now a constant quantity representing the distance from a constant phase 
plane to the screen. As shown in Figure (5), we make the projection of 
the field point r) on the screen, and denote it by O. Notice that O does 
not necessarily lie inside the aperture. Next, define two vectors I and If as 




(11) 




Field Point 



Figure (5) 
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shown in Figure (5). Then every point along I can be described by si, where 
< s < 1. Therefore 



f\ 



sT- If = ^S 2 l 2 + I 



where I and // represent the magnitude of I and If, respectively The area 
element on the aperture is 

da = Ids x sdl, 
and Equation ( TTTj) can be evaluated as 



1 

2^ 



aperture 



ik 



1 ^ e ikrl ±da 



2tt 



l f6 ik Ps 
2tt 



8=1 



ik 



edge J s=0 



(fx dlj 



s 2 l 2 + l) (s^P + lj) 



3/2 



e ik V s2l2+l f ffds x sdt\ ■ e z 



edge 



S 2 P + I) 



8=1 



S = 



But 



(fx df) 



where <fi is the angle subtended by the arc of the boundary of the aperture 
as measured from O. So 



ip(f f ) 



l f6 ik P s 

2tt 



Jkl f iky/P+Ff 



edge 



P + R 



1 

27 



i) s (f f ) -e ik ( ps+p f) cosOf 



where Of is the angle indicated in Figure (5). 

To separate the geometrical and reflected fields apart, we perform the line 
integral to the first term of the integrand: 



ip s (ff) d(f> 



ip s (ff) , if O lies inside the aperture 
, otherwise 



2" 71 " J edge 

Since O lies inside the aperture if and only if ff Gilluminated region, so 



ip(f f ) 



ip s (ff) , ff € illuminated region 
, otherwise 



2tt 



Jk(p s +p f ) 



COS 



edge 



(12) 
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Case 2. Point Source Diffraction 



Source Point 




Vj, Field Point 



Figure (6) 

Similar to the previous case, we attempt to assume that ff — f s is perpen- 
dicular to the screen, that is, we confine ff to lie on the central clXIS clS shown 
in Figure (6). But this time the assumption is totally unreasonable-after all 
we cannot restrict the position of ff. We believe that this formulation can be 
generalized without this assumption, but for the present we'll consider this 
special case, and try to elucidate the idea of "reflection at the boundary." 

Identify ip s = e tkp /p, where p = \\f— f s \\ , and define O, If, I as before. 
Here, we define a new vector l s , to be the vector from O to f s . Therefore 
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si — L 



s 2 P + l 2 s , 



and Equation. ( TTTj) can be evaluated as 



D ikp 



J aperture P 



ik 
r 



e lkr ^da 



8=1 



±f I 

2 71 " J edge Js=0 



ik( y /s 2 l 2 +l 2 f + y/s 2 l 2 +l 2 } 



ik 



27T Jedge 



Jky/sW+lj 



S 2 P + I) ( S 2 P + /2) 
e iky/s 2 l 2 +l 2 



3/2 



{ids x sdlj 



^2 + Z 2 J^2 + l 2 + ^/ S 2p + Ji 



s=l 



s=0 



So 



1 

2^ 



lf + ls 



Ps+Pf Pf 




e ifc(p s +p/) 

ib s (ff) cos 8f 
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where Of has the same definition as before. 

Again the geometrical field can be separated out by the same method, 
and we get the final result: 

, ( 4> s (ff) , ft e illuminated region 1 / e lk { Ps+p i) 

y t! I .otherwise ^J edgePs + p f iy 

(13) 



4 Solid- Angle Representation 
4.1 Motivation from Electrostatics 

Although Equation (JHJ) is a mathematically-consistent solution, it still ex- 
hibits the same problem as what Rubinowicz encountered in his solution: 
the geometrical and reflected fields are discontinuous at the boundary of the 
illuminated region. To overcome this problem, we seek for the analogy in 
electrostatics: Consider a grounded infinite conducting plane with a finite 
insulating region a at which the potential is held at a constant value Vq : 

y f Vo , on the insulating region a 
\ , on the conducting plane 

Assume there is no other source charge in the half space z > 0, and therefore 
V 2 V = there. This boundary value problem has the solution 

V (r» = ^V (14) 

where Q f is the solid angle subtended by the region o as observed at the field 
point ff. 

Inspired by the electrostatic result, we attempt a solution for diffraction 
problem of the form 

i\) (ft) = —^-ips iff) + (a line-integral), 

namely, we expect the geometrical field to take the similar form of Equation 
(|14|) . while the reflected field remains a line- integral around the boundary of 
the aperture. The advantage of this formulation is that both the geometrical 
and reflected fields now vary continuously, without any jump discontinuity 
across the boundary of the illuminated region. Different from Asvestas' work, 
in this formula, we see that as the field point ff approaches to the aperture, 
then £lf — > 2tt, and the geometrical field — > ip s while we expect the reflected 
field to vanish totally. That is, if we reside on the aperture, we are exposing 
ourself to the source ip s without the influence of the edge. 
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4.2 Inverse Cone 

We begin from Equation (J2J), and define 

J = -<f (^ s ^^-G K ^-\da, 

J aperture V @n On J 

J aperture V ulb un / 




Field Point 



Figure (7) 

To evaluate J, we do a trick slightly different from what Rubinowicz did. 
As shown in Figure (7), we make an auxiliary surface with the vertex at the 
field point, and make a small ball centered at the field point. Define o\ to be 
the surface of the cone outside the small ball while o"2 to be the surface of 
the small ball inside the cone. Apply the divergence theorem to the region 
enclosed by a±, <j<i and the aperture: 

F ■ da+ J F ■ da+ / F -da = 

aperture J a\ J cr'z 

where F = ipgVGx — Gk'V^'s, and da is the vectorial area element pointing 
outwardly from the volume enclosed. If we let the radius of the small ball 
approach to zero, then 

J.dd^-^(v.F)=-^ s (r f ), 
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where, as desired, f2 j is the solid angle subtended by the aperture as observed 
at the field point f F . So 



J = iff) - j F-da. 



<ji 



The surface integral can be evaluated by the same trick presented by Rubi- 
nowicz, as discuss in the following two cases: 

Case 1. Plane Wave Diffraction 

On the auxiliary surface, VGk = 0, and therefore 



r -, r p ikr 

/ F-da = / ike lkp p- 

J en J <ri T 



Pf x dl 



rdr 



Pf 



ike ik{r+p) Ps 



Pf x dl 
Pf 



dr 



ikp s 



p f x dl 



Pf 



--pf 



e ik{r+p) dr. 



From Figure (7), we have the relation 

P = ~(pf - r ) Ps ■ Pf + Ps 

and thus 



r r= Pf . . ., / „ . \ r=Pf / „ „ \ i J k (p s+ Pf) - e ik (-pfps-Pf+ps) 
/ e ik(r+ P )^ r _ e ik{-p fPs -pf+ Ps ) / e ( ps ' P! 'dr = " 

Jr=0 Jr= 



'r=0 



ik 



and 

F-da 
Therefore 



I F-dd=l e* kp ° e — (l - e-^(^ p f)) ( 8l2L£L.) 

Jen /edge Pf V ' V 1 +Ps 'Pf) 



1 + Ps-Pf 



dl. 



J 



— fa (r f ) * 



An 



Jkp s 



edge 



Jkp f 
Pf 



-ikpf(l+p s -Pf) 



Ps X Pf 
1 + Ps-Pf 



dl. 



To evaluate J*, we use the result from Section 2: 



J 4tt£, 



D ikp s 



-,ikpf 



Ps X p* f 



'edge pf yt + Ps-p} 

Finally, we combine J and J*: 



■dl. 



i> iff) = ^s {ff) 



4vr/ cdgc 



Jkp s 



Jkpf 
Pf 



, _ e -ikp f (i+p s .p f )\ ( Ps x Pf 

J V 1 + ps ■ Pf 



Ps X p) 
1+Ps-P*f 

(15) 



dl. 
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Although this result fits our demand-?/' (rf) is now expressed in terms 
of the solid angle fiy-it is still unsatisfactory since the denominator of the 
geometrical part is 4tt instead of 2tt. Accordingly, if r/ approaches to the 
aperture, the geometrical part only gives us one half of the total source wave 
ip s , and thus the reflected part must contribute the rest half part. To fix the 
problem, we take the long wavelength limit k — > 0, and thus 

^ e ^op = 1; 
and Equation (fl5j) must be identical to Equation (fT4l) : 

Ws)^ + H -( psXp \)-diA 

VKIJ 4vr ATrJ edgePf \l + p s -p* f ) 2tt 
So we have a line-integral representation of solid angle: 

^ = ±t \.£ (i6) 

This equation has also been derived by Asvestas [3] , and by Yih-Yuh Chen [I] 
from a more elegant perspective. Note that since the wave number k vanishes, 

the vector p s now can point in an arbitrary direction, so the representation 
above is not unique. 

Finally, we construct the desired ip iff) by adding Equation (1TB]) into 
Equation (|T5|) 



Vtt 1 / e ikpf 

2vr 47ry edge p f 



Case 2. Point Source Diffraction 

Again, on the auxiliary surface, VGk = 0, and 




■all. 



F.da= I t —^(^-^\p.[ P -t^rdr\= / e**+» 




n 2 



p PV V Pf J J*i \P p 



Ps e ik(r+p) (te_}_\ dr 



r=0 \P 2 P 3 
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Field Point 



Figure (8) 
From Figure (8), we have the relation 

P 2 = Pl + (P/ ~ rf - 2p s (p f -r)p s - p f . 
Differentiate it with respect to r : 

dp 



P ' 1+ dr~> =r + P-Pf + P*P^Pfi 



and (follow Rubinowicz' calculation) 



,ik(r+p) 



r=0 



ik 1 

p 2 p 3 



dr 



3 ifc(r+p) 



p{r + p-p f + p s p s ■ pf) 



=/-'/ 



r=0 



Jk(p f +p a ) 



pl (i + Ps ■ Pf) pI (i + po ■ Pf) 



where p = p s — pj is the vector from r s to ff. So 

- -|— ' 

edge 



F ■ da 

I (7 1 

Therefore 

n 



2 „-ik(p s +pf-po) 
rs c 



Ps Pf [l+Ps'Pf Po 1 + PO-Pf 



(Ps X Pf) ■ dl 



J = gV- (ff) 



-Li 



Jk Ps e ikp f I j p2 e -ik(ps+pf-po) 

47T 7 cdgc p s pf \ 1 + Ps ■ Pf Po 1 + PO ' Af 

Again, we use the result from Section 2: 



(p 8 x py) • dl 



J* 



47T 



; ifcp s e ikp f / p s X p* f 



Ps Pf ^1 + pVP/ 
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Combine J and J*, we have 



ij;(f f ) 



5/ 

47T 



■y r e ikp s e ikpf 
47r 7cdgc Ps Pf 



1 



2 -ik(p s +pf-po) 



l + Ps'Pf Po 1 + PO ■ P/ 
Ps X P*f 



(Ps X Pf) 



1 + P s • P* f 



■dl 



(17) 



To construct the correct factor 1/27T, we use Equation ( TT6|) to add another 
(n//47r) ^ s to the geometric wave. But note that p s in Equation (TIB]) is 
an arbitrary constant vector, and in Equation (11 7p p s represents a varying 
vector that changes its direction as we integrate along the edge. Thus we 
must specify one direction for p s in Equation ffT6]) so that we can insert it 
into Equation ffTT|) . The result is most symmetric if we adopt p s = p in 
Equation (fTB]) 



ikpf 



2tt 47T7 edge p f 



2 p -ik(p s +Pf-po) 
rs 



1 + Ps ■ Pf Po 1 + PO • Pf 



(Ps X 



Ps X Pf 
1 + Ps ■ P*f 



_|_ e -ikp f 



PO X p / 

1 + Po ■ p} 



We have seen that in both cases the field ip (ff) has the form 



—ips (r/J * 



2vr 



edge 



P.f 



and it is obvious that both the geometrical and the line integral parts of ij) 
are now continuous across the surface of illuminated region. However, the 
integrand in [• • • ] now depends on the type of the source. There is another 
point to be mentioned: since Qf oc 1/p/ 2 , in the far zone the geometrical 
field is overwhelmed by the reflected field, which is proportional to 1/p/ of 
the source wave. The situation is reversed in the near zone, of course. 



5 Comparison of Boundary Conditions 

As discussed in previous sections, the boundary conditions based on Kirch- 
hoff 's theory is mathematically inconsistent, and by using the proper Green's 
function, the diffraction theory can be transformed into a boundary value 
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problem of Dirichlet type which is mathematically admissible. However, this 
is not the whole story. As the source wave ip s propagates toward the aper- 
ture, the wave must be modified by the presence of the opaque screen, and 
thus if) is not exactly equal to if> s , the unperturbed source, on the aperture. 
So the boundary values Equation (J4]) imposed earlier is still, unsatisfactory 
in the physical sense. 

However, Sommerfeld has solved a 2 — D straight edge diffraction problem 
rigorously without using the unperturbed source wave as boundary values [5], 
and we'll see in a moment that, by some deformation of Equation ([8]), the 
functional form of our solution is very closed to that of Sommerfeld's, and, 
therefore, we may regard Equation (j4]) as an acceptable approximation to 
the real, rigorous solution. 



5.1 Approximated Solution for a Point Source 



Consider an infinite half plane lying on z = and x > 0, with a point source 
lying in the region z < as before. The solution if> in the space z > can 
be solved by Equation (jSJ) 



if;(f f ) 



if> s (ff) , ff e illuminated region 
, otherwise 



1 



A 



e lk Ps e lkp; I p g X f)f 



Ps X p) 



■ dl. 



4^ ./edge Ps Pf \ 1 + Ps ■ Pf 1 + Ps ■ Pf 

where the line integral is performed along the infinite straight edgd3- The 
amplitude A of the source field is now expressed explicitly for later conve- 
nience. 

In the far field region kpf ^> 1, we apply stationary-phase approximation 
to evaluate the reflected field: 



An 



ik(ps+Pf) 

A- 

edge PsPf 



Ps X Pf 



Ps X p f 



l + p s -pf l + p s -pf 



■ dl. 



;is) 



The the stationary point occurs when V (p s + pf) • dl — 0, and we expand 
the phase at the stationary point: 



Ps + Pf = {Ps + Pf) 

= (Ps + Pf) 



+ V (p s + Pf) 



51 + -61 



PsPs ]_ 



pfpf 



pf 



51 



, 1 Ps + Pf r?l . 2 1 ,, 

H -51 sin [dl,p< 

o 2 PsPf 



o 

(19) 



1 And thus the integral sign J is used instead of §. 
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where the subscript denotes the stationary point, which in this special case 



is the point on the edge nearest to f/; also, 1 is the identity operator in three 
dimensional space, and sin (dl, p s j is the sine of the angle between dl and p s . 
Insert Equation (TT9~j) into Equation (IT5|) . and perform the Gaussian integral, 
we get 



1 Ae ik y ps+Pf ' 



Ps X pf Ps X p f 



dl 



4vr p s p f \l + p s -pf l + p s -p* f 



dl 



2ttz p s p f 



1 



k P' + Pfwn(cU,p 



But 



Ps X Pf Ps* Pf 



dl 



1 + Pa-pf l+p s -p* f 



dl 



p r n 



Pf 



n 



l+p a -pf l + p s -p* f 



- sin [dl,p 



where 



n 



dl x p s 



dl x p. 



is the unit outward normal of the geometric light cone. To simplify the factor 
in the parenthesis, refer to Figure (9), we have 



p r n 



Pf' 



sin 



a 



sin [6 + a) 



2 sin i 



1 + p s • Pf 1 + Ps • Pf 1 + cos (0 — a) 1 + cos ((f) + a) cos a + cos i 




Figure (9) 
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So finally, 



-A 



Jk(p s +p f ) 



2 sin< 



e * 



cos a + cos y/pT+~Pf 



(20) 



The factor e*4 can explain the reason why the diffraction pattern in the water 
has a phase delay compared to the incident wave. 

5.2 Approximated Solution for Plane Waves 

In the case of plane- wave- incidence, we can obtain the solution from Equation 
(120]) by taking the limit 

p s — > oo, A — > oo, while keeping > finite number, taken to be 1 

Ps 

The result is 



Jk(p a +p f ) 



2 sin< 



^/2-Kkpf \ cosa + COS( 



e *i = L±L e ^ 



cos • 



cos ■ 



COS 



sin 



Here e* fcps has been dropped since Sommerfeld assumed that the plane 
wave has phase right at pf — 0. 
So the total field is 




4> s (f f ) + l +^ e ik Pf 
W nk Pf 
l + i 



In comparison with Sommerfeld's solution [5] 



cos f 

— , ff G illuminated region 

sin | 



cos f 

otherwise 



sin 



. . l + i ., 
^ (ff) + -^=e lk Pf 

W nk Pf 

l + i 
4^/nkpf 




, ff G illuminated region 
, otherwise 



we see that, apart from the factor 



cos f 
sm | 



the two representations are similar. 



The discrepancy results from the different boundary conditions as discussed 
before. 
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If, however, we use the Kirchhoff 's integral formula (with mathematically 
inconsistent B.C.s), we would obtain 



which has a different functional form from Sommerfeld's solution. 

6 Conclusion 

By using mathematically consistent boundary conditions, we have seen that 
Rubinowicz' decomposition formulation can be more useful: the functional 
form of the line integral becomes much neater and admits a simple inter- 
pretation of reflection at edges. The formulation also provides us a different 
approach that makes the diffraction phenomena similar to the electrostatic 
problem by using solid angle representation for the geometrical field. Fi- 
nally, the diffracted field predicted by this formulation is much closer to the 
physical solution, as discussed in the last section. 
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